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ABSTRACT 

We study the properties of density perturbations of a two-component plasma with a tempera- 
ture difference on a homogeneous and isotropic background. For this purpose we extend the 
general relativistic gauge invariant covariant (Gic) perturbation theory to include a multi-fluid 
with a particular equations of state (ideal gas) and imperfect fluid terms due to the relative 
energy flux between the two species. We derive closed sets of GIC vector and subsequently 
scalar evolution equations. We then investigate solutions in different regimes of interest. In 
particular, we study long wavelength and arbitrary wavelength Langmuir and ion-acoustic 
perturbations. The harmonic oscillations are superposed on a Jeans type instability. We find 
a generalised Jeans criterion for collapse in a two-temperature plasma, which states that the 
species with the largest sound velocity determines the Jeans wavelength. Furthermore, we find 
that within the limit for gravitational collapse, initial perturbations in either the total density or 
charge density lead to a growth in the initial temperature difference. These results are relevant 
for the basic understanding of the evolution of inhomogeneities in cosmological models. 

Key words: cosmology: theory - plasmas - relativity - gravitation - cosmology: large-scale 
structure of Universe 



1 INTRODUCTION 

Plasmas and electromagnetic fields are common in our Universe. 
They play an important role in a diverse setting of astrophysical 
and cosmological processes. Plasmas may be found, e.g. in stars, 
accretion disks of rotating black holes, the Earth's ionosphere, and 
also constitute the intergalactic, interstellar as well as the intrasolar 
medium. There also are occasions when general relativistic gravity 
has to be taken into account in conjunction with plasma physics, 
such as in the close vicinity of the aforementioned rotating black 
holes. Another prominent example is our Universe, in which the 
plasma state has over time been more or less prominent. Obviously, 
in such situations gravitational effects due to general relativity, such 
as gravitational waves, can not be neglected, and may lead to inter- 
esting results. 

Moreover, perturbation theory within gravitational physics has 
always held a special place within physics, due to, e.g. the pioneer- 
ing work of Jeans, and its relation to the basic question of structure 
formation. The issue of structure formation is now a mature science 
and cosmological observations have been taken into the high preci- 
sion regime by observational tools such as COBE, WMAP, Chandra, 
and the planned LISA mission The development of new observa- 
tional instruments also calls for new theoretical models and con- 
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cepts to be tested, in order to refine the current standard model of 
cosmology. Thus, there is ample interest in extending the existing 
models by including other physical effects. 

In perturbation theory there are two main distinct schools 
of thought. The first could be described as metric based and fol- 
lows from a seminal paper by Bardeen ( 1980) which was extended 
to multi-fluids bv lKodama & Sasaki ( 1984) a nd dev eloper further 
in a more recent paper by .Malik & Wanda i2005h . The second 
paradi gm is a covarian t an d s^auj^e -invariant pe rturbation theory 
due to Hawkine f 196d) and Ellis & Bruni fl98?) (for reviews see 
Ellis 1 1995) or EOis & van Elst 1 1 999)), which was extendgd t o 
multifluids in Dunsbv et al. ( 199J): lMarklund et alj J200(l 12003) : 
■Betschart et al.. ^2004,) and applied to the cosm ic microwave back- 
ground (CMB) bv lChallinor & La senbyHl 998h. The equ ivalence of 
the two theories has been shown in lPunsbv et alJil992h . 

We will adopt the second approach because the physical quan- 
tities defined as gauge invariant and covariant (GIC) variables are 
not only convenient mathematically, because they avoid spurious 
non-physical gauge modes, but also clearly reflect the physical 
quantities that an observer would measure and these variables al- 
low a clear geometrical interpretation. In this paper we will apply 
the theory to a two-component, two-temperature plasma in a cos- 
mological setting. In particular, we assume that the two plasmas 
each satisfy the perfect gas law, but with different temperatures. We 
formulate the set of governing equations for the plasma dynamics 
on a general relativistic background. These equations are analysed 
by perturbing the two-fluid model around an isotropic and homo- 
geneous background using the covariant gauge invariant approach. 
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A large-scale temperature difference between the electron and 
proton fluids in a cosmological setting is not very likely since the 
Coulomb equilibration time is usually much shorter than the growth 
rate of perturbations. However, local fluctuations might occur and 
more importantly we hope to apply the generalized two-fluid equa- 
tions derived in this paper to other two-temperature fluids in fu- 
ture papers, such as dark matter - ordinary matter, decoupled free- 
streaming neutrinos in the early universe and possibly the advection 
dominated accretion flows onto black holes, where the electrons 
and protons are thought to be out of thermal equilibrium. 



2 GAUGE INVARIANT COVARIANT THEORY 

To make this paper more self-consistent we will recapitulate some 
of the results in the literature that we use here and make the appro- 
priate approximations to linear perturbations. 



2.1 3+1 split 

To start with, we will consider two perfect fluid species that flow 
in a curved space-time. Both in an astrophysical and a cosmologi- 
cal context it is often possible to identify a prefeiTed family (con- 
gruence) of fundamental world lines associated with the motion of 
typical observers. In the former case, one could for example choose 
the world lines of observers that move with the bulk velocity of the 
matter in a (relativistic) plasma wind or jet associated with a pulsar 
or active galactic nucleus. In cosmology it is customary to treat the 
Universe as consisting of a perfect fluid whose world lines are de- 
termined by the motion of distance clusters of galaxies with respect 
to an observer. 

When such fundamental word lines can be identified, it is pos- 
sible and convenient to split 4-dimensional space-time into 'time' 
and 'space' again with respect to the 4-velocity vector tangent to 
the world lines 



dx" 
1x 



-1 



(1) 



where T is the proper time along the world lines. The tensors that 
project onto u" and into the tangent 3-spaces orthogonal to u" are 
given, respectively, by 



U"), = -«"«/, and hat, = gab + «fl"i 



(2) 



Associated with the 3-1-1 split are two derivatives. The 'time deriva- 
tive' is the covariant derivative along the fundamental world-lines 
and the 'spatial derivative' is the covariant derivative of any tensor, 
where all the indices are projected on the hypersurface orthogonal 
to u". For an arbitrary tensor T"*^ these are defined by 



hab 



'V.T 



ab 



cd 



v„r 



ab 



cd 



= hyi%hP^h%h'-^W,T 



/S 



(3a) 
(3b) 



The derivative V is a proper 3-dimensional derivative if and only 
if the vorticity of i/' is zero, for instance in a Friedman-Lemaitre- 
Robertson-Walker (FLRW) Universe (or a FRW Universe when we 
neglect the cosmological constant). 



2.2 Multi-fluid stress-energy tensor 

We assume a congruence if corresponding to fundamental ob- 
servers as defined in Section lOl and a Universe filled with a plasma 



made up of two species that are allowed to move in arbitrary direc- 
tions with respect to the observer with individual (non-relativistic) 
4-velocities satisfying 



"(0 



-W+v- 



(4) 



The index ; = 1, 2 or ; = +, — labels the two fluids. Each species 
is assumed to be a perfect fluid in its own rest frame. Then m?-, is 
the unique hydrodynamical 4-velocity that is time-like and allows 
a split of the energy-momentum tensor in the perfect fluid form' 



rab 



^M(,: 



-P{' 



,A' 



,ah 



(5) 



Here p/,) is the pressure and /i/,^ the energy density and for a perfect 
fluid M?., is also parallel to the particle and entropy flux, such that 



Nf:, 



•i(,-)!, 



Si: 



'■)■ 



\i) 



(6) 



and there is neither a particle drift nor an energy flux (J?.. = cfi.s = 0, 
nif. the number density and 5(,) the entropy). 

In the frame of the fundamental observers ;/", which in gen- 
eral will move relative to the rest frame of the individual fluids, 
the energy-momentum tensor is clearly different from J5j. In par- 
ticular it can be written in the general form of an imperfect fluid. 
Eq. ^4} defines the velocity v?., of the fluid i with respect to the ob- 
server u" . If we assume that in the background all species share the 
same hydrodynamical 4-velocity m" then v?., = in the background 
and is therefore gauge invariant by the lemma of lStewart & Walken 
( 1974). 

We linearize the stress-energy tensor for non-relativistic devi- 
ations from the average flow such that \^,.. <^ 1 and the an-isotropic 
pressure is negligible, but to linear order there is a contribution to 
the energy flux, or heat, given by cfi.s = {piu) +P(,))v?> such that 



r(^,^ = M(o""«'' + P(,y^'" + 2"'''4' 



(7) 



which are the linearized equivalents of Eqs. (3-4) in lMarklund et alJ 
<2003h. 



2.3 Equation of state 

In this paper we will use a particular equation of state to study 
temperature effects in a two-component plasma by a first order per- 
turbation analysis on a FRW background. We assume a ideal gas 
equation of state for each species and corresponding energy den- 
sity and enthalpy hij\ : 



P(i) 
M(,-) 

h(:\ 



(™(i) + i^«^(o)"(0 ' 

^{i)+P(i) = (™(f) + i^B^(/)) «(/) ' 



(8a) 
(8b) 
(8c) 



where Tii\ is the temperature, mn\ the rest mass, and kg Boltz- 
mann's constant. We will allow for different temperatures between 
the species, and note that Tn\ is time-dependent in the background 
FRW space-time. Moreover, we will re-scale ksTn) -^ Tn) for the 
sake of brevity. 



This is the u| 4-velocity in lPunsbv et alj 1 19921) . 
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2.4 Conservation of particles, energy and momentum 

Conservation of energy and momentum is guaranteed automati- 
cally by the twice-contracted Bianchi identities and Einstein's field 
equations and is given for each fluid by 



v,r«)* = «,F«4, 



(9) 



where T,".!' is the stress-energy tensor for the matter and the right- 
hand-side follows from the electromagnetic part by V^T — 



em.(/) 



—IM)JP%J() with 7?.> the 4-current density and /Jo the vacuum mag- 
netic permeability. Similarly, conservation of number density of 
particles is given by VaA'fx = 0. 

The time-like and space-like projections of J9l provide 
separate linearized equations for e nergy and momentum conserva- 
tion, respectively (see for instance lMarklund et alJi200Cl) ) 

A(,-) = -%)(0 + V„v^,.)), (10a) 

/;(,■) (m^ + vJJ) = -VV(,-) - (f /^(,-) +P(/))^(/) +P{i)E" , (10b) 

where is the expansion, pj,-) = q(i)nij\ is the charge density and 
E" the electric field strength. 

A separate evolution equation for the temperature can be 
found from combining <8bt and JlOat with the linearized equation 
of particle conservation: 

n(,-) = -(0 + V«v;',.) )«(,.), (11a) 

^('■) = -f(® + Mo)^«- (lib) 

2.5 Expansion 

The basic equation for gravitational attraction is the Raychaudhuri 
equation that in a FRW Universe with vanishing acceleration, shear, 
vorticity and cosmological constant^ reduces to 






(12) 



with K = 87i:G/c'*. lEiiis & van Els3 h999tl call the term /i + 3p in 
<12t the active gravitational mass density of the plasma to empha- 
sise that this term and not just fl drives gravitational contraction. 
The volume rate of expansion = Vhk" is the trace of the covari- 
ant derivative of u" and is related to the Hubble parameter H and 
the scale-factor S by 



Sit) 
S{t) 







H . 



(13) 



2.6 Electromagnetic field equations 

The electromagnetic field is determined by the linearized Maxwell 
evolution and constraint equations: 



£("> = -lOE" + curlS" - flof 



V„B" = , 



|0S°-curl£"' 



(14a) 
(14b) 
(14c) 
(14d) 



All equations can easily be generalised to include a cosmological con- 
stant, but for simplicity we will specify to a FRW Universe. 



Table 1. Variables for the summed and subtracted quantities for the two 
fluid species. 







Total 


Differential 


Energy dens. 


M 


= M(1)+M(2) 


5ai = M(i)-M(2) 


Pressure 


P- 


= P(1)+P(2) 


Sp = P{l)-P{2) 


Enthalpy 


h 


= /l(l)+ft(2) 


5/! = /l(,)-/l(2) 


Temperature 


T 


= T(l) + Tf2) 


ST = 1(1) - 7(2, 


Number dens. 


N-- 


= W(l)+%) 


5" = "(l)-«(2) 


Velocity 


V = 


5(f(l)+^'(2)) 


5v=5(v(i)-V(2)) 


Mass 


M-- 


= m(i)+m(2) 


dm — mn\ —ffln) 



where p = g(i)n(i) +?(2)n(2) and / = ^(ijnjijv^jj +<?(2)"(2)Vp) 
are the total charge density and total 3-current density, respectively, 
qn\ is the particle charge and Eq the vacuum electric permittivity. 

To be more specific, we will choose species 1 to correspond to 
protons or positrons with qn■^ = e and species 2 to electrons ^o) = 



2.7 Definition of gauge invariant quantities 

Since the background space-time is homogeneous and isotropic, we 
require that the electromagnetic fields vanish to zeroth order This 
implies that the charge density and current density should vanish 
and thus p=0=>n(i)=n(2)' and j" = enM^ . — v^^ ) = => v?j > ^ 

^(2) 



Because the energy flux should also vanish in the background 



^^2) 



= 0. 



Thus, the only nonzero variables in the background are the ex- 
pansion 0, the energy density /x,,), the number density M,,, and the 
pressure pij\, or equivalently the temperature Tuy The background 
evolution of these quantities follows from the zeroth order terms in 
^, (TOa} and (TT} 



f^(i)=-®{P-{i)+P{i)) > 



^(0 
Pii) 



=-i0p(/) 



(15a) 
(15b) 
(15c) 



The remaining variables are gauge invariant covariant (GIC) 
first or der perturbations on the FRW background by the lemma in 
IStewart & Walker 1 1974) . For an extensive physical motivation of 
using gaug e invariant variables that vanish in the (fictitious) back- 
ground see lBardeerJ h980l) : lEllis & Brunil J1989h : lEllis & van ElstI 
( 1999). 

In the next section we will define GIC variables for all the 
physical quantities that govern the two fluids separately, and their 
mutual interactions and derive covariant evolution equations for all 
these variables. 



3 EVOLUTION OF GIC VECTOR PERTURBATIONS 

We now proceed to find evolution equations for the combined two- 
fluid as determined by the observers u" (note that until now this 
frame is still left arbitrary). To facilitate this we define the variables 
in Table[T] 

In terms of the variables in TableQwe can write the first order 
charge and current densities as: 



eSn . 



eNdv" 



(16) 
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Or alternatively, the velocity difference can be interpreted as a nor- 
malised current density dV = j''/{eN) and similarly, the number 
density difference Sn is clearly equivalent to the charge density 
8n = p/e. 

The evolution equations for the density and temperature vari- 
ables (with 5 r 7^ in the background) follow from Jilt and <15bt : 



Table 2. GIC dimensionless comoving scalar variables. 



A' A' 

A' 

St t - ST T ~ 



(17a) 
(17b) 
(17c) 
(17d) 



Similarly, we find the total momentum conservation equation 
by summing <10b> over the two species 

fe(" = -V°p-|0^"-4<°> , (18) 

while subtracting <10bt yields 

5hv'''''>+h5v''"^ = -{§Sh + 8p)V" - (§h + p)Sv" 

-8hu"-V"Sp-5pE" . (19) 



3.1 Frame choice 

From hereon it will be convenient to specify to a particular ob- 
server frame u" . Customary choices are either the particle frame in 
which to linear order L"(i)^'(i) = ^^" = such that the total ve- 
locity perturbation with respect to the observer always vanishes, or 
the energy frame in which q" = Y,h(i)'^(i) = ^^" + ShSv" = 0. 

Both are useful in eliminating either V" or Sv". In the par- 
ticle frame, this relation is particularly simple with V" = and 
5v = V(i) = V(2)- However, requiring the total velocity perturba- 
tion to vanish when studying wavelength dependent perturbations 
in the total density, for instance, doesn't seem like an obvious frame 
choice. Choosing the energy frame means that the acceleration is 
driven only by pressure gradients: hu" = —V"p from the total mo- 
mentum equation <18> . The relation V" = dv"{Sh/h) is convenient 
in a cold plasma when /ij,) = miisni^-s and 8h/h is time independent. 
In our case, however, when the pressure contributes to the energy 
density, this coefficient has a complicated temporal dependence and 
doesn't really simplify the calculations. 

In this paper we therefore prefer to choose the geodesic frame 
of freely-falling observers such that if = (and cf ^ 0). The cor- 
responding 4-velocity is covariantly defined and reduces to the per- 
fect fluid flow velocity of the fundamental observers in the back- 
ground. Therefore the GI variables defined by gradients orthogonal 
to the fluid flow are properly defined and the geodesic frame is a 
valid choice. 

One could interpret <18t as defining Sv" as a function of V" 
and eliminate V in favour of Sv" in <19t to get an equation similar 
to those in the energy frame. Rather, we treat <18> and <19t on equal 
footing as two coupled differential equation for the total velocity 
and the normalised current density. By eliminating p and 8p using 
<15c> . we find from <18> and <19> the (almost) symmetrical relations 



Total 


Differential 


X" = S^ 


y=s^ 


y-s^"/ 


f-S^°f 


W' = SV^VfeV* 


vy° = 5V''V45v* 


Z" = S^"@ 





/;V'<"> + dhdv'^"^ = -®{^-Ap)V" -§{Sfl-45p)Sv" 

- V"p . (20a) 

5hv''"'> +h5v^''^ = -§{5fl-45p)V" - §{fl-4p)dv'' 

-V^Sp-dpE". (20b) 

The above equations are clearly GIC since V" ,Sv" .^"p, V''Sp and 
E" all vanish in the background (implying that all the coefficients 
in <20t should be evaluated in the background). 



3.2 Definition of GIC vector quantities 

It is argued in lEllis & BruruHl989l) that the physically relevant vari- 
ables characterising the spatial variation of scalar quantities that do 
not vanish in the background (such as the energy density) are the 
dimensionless comoving fractional spatial gradients of those quan- 
tities. In particular, we define GIC variables for the density, temper- 
ature, expansion and velocity perturbations in Table|2| 

The variables such as X" are the quantities that a typical ob- 
server would measure. For instance, at a particular time one can 
measure the comoving gradient in the number density (with re- 
spect to the length-scale given by the background scale-factor S) 
by simply counting the number of distant galaxies or clusters of 
galaxies at each distance (obviously, closer by sources like those 
in our own galaxy do not follow the 'cosmic fluid' and are neither 
homogeneous nor isotropic). Other choices like 5/x//x comparing 
the perturbed energy density with its background value depend on 
the fictitious background and are gauge dependent to the extent that 
some perturbations can be eliminated by a choice of gauge and the 
distinction between physical modes and spurious gauge modes can 
be obscured. 

The velocity variables such as V" (V;,l/ ) might be interpreted 
as the gradient of some kind of fluid compression, but in this paper 
will be treated simply as auxiliary variables needed to couple J20> 
to(T7}. 

3.3 Evolution of GIC vector perturbations 

To obtain GIC evolution equations for vector perturbations in the 
total and relative densities and temperatures and the expansion we 
take the comoving gradients of <17t which conveniently commute 
with the time derivative (this follows from the commutation rela- 
tions in the Appendix <AH and the definition J13» . The expansion 
terms proportional to in <17> are absorbed by choosing dimen- 
sionless variables. For example from <17at and <17bt we find (to 
linear order): 

Sn \ Sn „ Sn 

N J N N 

The resulting GIC evolution equations for the vector quantities 
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in Table|2|are 

(21a) 
(21b) 
(21c) 
(21d) 

(21e) 

where m = 5m/M, T = T/M and a = ST/T. Note that ^ are 
still valid in any observer frame u" and also that a is constant in 
the background. 

We define Q." = curlV, a" = curl5v°, and "^^ = curl£° and 
take the curl of <20t . The pressure gradients vanish and using <AU 
we find 



Table 3. Gauge invariant scalar expansion normalized variables. 



x<«> 


= 


-[W"+Z"] , 


^{0) 


= 


-w" , 


y(«> 


= 


-l[W" + aw"+Z"] , 


yi") 


= 


-^[aW" + w'' + aZ"] , 


2{a) 


= 


20Z« MN ^9fY. 




+ 


(l + §f)X'' + (m+|af)x"; 



hCl" + Sha" = & [(p- |M)fl" + {dp- |5/i)ffl"] , 

hco" + ShQ." = @[{p- |m)g)" + {Sp- |5/i)fl"] +eN'^" 



(22) 



The curl-free part of the velocity perturbations follows from 
taking the divergence of <20t and eliminating VaE" in favour of a 
charge density perturbation and the background plasma frequency 



"P-O) 



e 



"(1) I "(2) 



^(1) 



^(2) 



eoMl- 



(23) 



More important though is the gradient of the resulting equations, 
which yields the equations for W and w" that couple back to the 
density and temperature perturbations <21> . These GIC auxiliary ve- 
locity variables evolve as: 

hw'^"^ + Shw'^"^ = 

- ^ [v2 + I (©2 - 3m)] [X" + Y" + ax"] (24a) 

-&[{lti-p)W" + {l8ii-8p)w"] , 

and 

-^[v2 + |(02-3/x)][aX«+/+x"l 
-@[{y-p)w" + {p^i-Sp)W"] , (24b) 

where we have used J14ct . the commutation relations JA1> and 



NT 

-^{X" + ax"+Y") 



V"8p = —{aX"+/'+y^) 



(25a) 



(25b) 



3.4 GI electromagnetic equations 

To investigate how the velocity perturbations couple to the elec- 
tromagnetic field we need the solenoidal part of <2Qb> . which cou- 
ples to B" through <14b> and the curl-free scalar part which through 
<14ct couples to the plasma frequency <23t . 

From the constraint <14ct and <21bt we trivially find an evo- 
lution equation related to the curl free part of the electric field in 
terms of SV"(Vi£'') = ^° = eNx"/Eo 



<r<°>: 



eN 

> 

£0 



(26) 



Total 


Differential 


A = 52^ 


s-s'-%'" 


^=s^^L 


r=s'-^ 


#- = S^V^VfcV* 


r = s^v^Vi,5v'' 


f =S-V20 





y2Qa)) the wave equations to linear order 

^{a)_y2^a ^ I (^ + 3p)£"-f ©£:<''> -i02£« 

- — V^Sn- floeNSv'"''^ , (27a) 

£() 

g(fl>_v2gfl ^ ^(^ + 3p)S°-§eB<">-^0V 

+^eNco" . (27b) 

The expansion normalised curl of the electric field '/f" = S curlE" 
satisfies the wave equation 

^■(«> _ vV" + fO'^^''^ = i(/i -f 3p - |02)^° 

-^eA'[co<°> + fffl''], (28) 

from i27al and JaT}. 

The evolution equation derived in this section together with 
the constraints completely determine the behaviour of GIC vector 
perturbations of the two-temperature plasma. Here we will focus 
on scalar perturbations. In the next section we will derive the cor- 
responding evolution equations. 



From Maxwell's equations, we find (in agreement with iTsagaj 



4 GIC SCALAR PERTURBATIONS 

The easiest way to obtain a closed set of differential equations for 
scalar perturbations of the two-temperature plasma is to take the 
divergence of the vector equations derived in the previous section. 
The resulting GIC expansion normalised dimensionless variables 
are summarised in Table|3| 

In section IT2I we argued that the variables in Table |2| which 
are dimensionless fractional spatial gradients of scalars quantities 
that do not vanish in the background space-time, are in fact the 
physically relevant variables that an observer would measure di- 
rectly. By taking another divergence of these variables and subse- 
quently performing a harmonic decomposition (in the next section), 
we essentially introduce a length-scale V ~ L^ corresponding to 
a Fourier component of a perturbation with wavenumber k = 2n/L. 
Thus A, for instance, represents a comoving density perturbation on 
a length-scale L: A --^ iS/L) [dp / p), which is similar to the vari- 
able 8p/p that is sometimes used, but defined to be covariant and 
gauge invariant. A similar interpretation applies to the other vari- 
ables in Table|3| 

Since the comoving divergence again commutes with the time 
derivative, the scalar equations look very similar to <21L All the 
coefficients have to be zeroth order and it is more illustrative to re- 
write those in terms of the background quantities A', T, p,„ = MN /2 
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etc. The resulting evolution equations are 

A = -[#- + ?], 
8 = -r , 

f = -^[W + a^r + Q = l[a8 + A], 



(29a) 
(29b) 
(29c) 
(29d) 



(l + §f)A- 



rr 

9„T■^ 



;+faf)5] . (29e) 

In com paring with the non-thermal dust results in lSetschatt et alJ 
J2004 we find that i29bll29al agree with their (24a,d). For the lat- 
ter identification one should transform back from the energy frame 
where V" = —Sv^dfi/fi was used to eliminate V, and instead 
leave the equation in terms of 5v". Clearly the equations for the 
thermal perturbations are new and the evolution equations for the 
velocity and expansion perturbations are more complicated when 
including temperature effects. 



4.1 Harmonic decomposition 

In terms of the variables in Table|3| the evolution of the scalar ve- 
locity variables is found by taking the divergence of <24all24b1 us- 
ing lAH .-^ Furthermore, we harmonically decompose the Lapla- 
cians in terms of the scalar harmonic functions g' ' that satisfy the 
Helmholtz equation 



y2Q{i) . 



'52 



Q 



W 



e(^'' = o, / = E/(,)e<' 



(30) 



where f is any perturb ation scalar and f,/,-, the corresponding har- 
monic iHarrisonlll967l) . For the sake of clarity we will drop the 
index £, as we will not encounter harmonic mixing due to the lin- 
earity of the equations. Each evolution equation is understood to 
represent the temporal behaviour in one harmonic mode £. The co- 
nioving wavelength of a perturbation is given in terms of the wave- 
number £ by A(t) = 27tS{t)/£ where time dependencies are given 
in terms of T = t/tQ, the proper time along a world-line normalised 
to some initial time fg- The velocity variables then evolve as 



l + ff)#+(m + faf)f : 



(31a) 



£^T 



^{A + aS+T)-^@{W + mf) 



l + ^f) f+(jn+^aT]W = 

£^T , 



(31b) 



'-^{aA+8 + y)-^{r + mW) + {l-,fp-)cof,8 



4.2 Sound velocities 

We assume that the background plasma is non-relativistic such that 
the total thermal energy is small compared to the total rest-mass 
energy. Consequently, T = T /M <ti 1 and we will neglect f in 
the gravitational interactions and use /x ~ ;U + p ~ Pm in the back- 
ground. 

We define the sound velocity for each species as: 



a--) 



A(/) 



-n{i]T[i) 



57» 
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(32) 



The commutation rule for V„V- cancels the factor proportional to the 
expansion in the previous V"V^ commutator. 



using <15> . We define total and differential sound speeds similar to 
V" and Sv" as 



C 



5c? 



(cL+42) 



jy-sA 
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I -ma 
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:±ff 



(33a) 



(33b) 



where the sign of <33b> is positive when c^ j > c^2 ^^'^ negative 
when c^ J < c^j. The variable Cj is related to the physical mass- 
weighted total sound velocity c^ = Y,(i) c,..iin)/ii = 5p/3fl (for our 
particular equation of state) by 

r2 _ 2 f l-{Sfi/fi){5p/p) \ _ 2 f l-ma \ 

where the second expression (Cj = 5r/3) follows from the fact that 
the number densities are the same in the background. In the remain- 
der of this paper it will be apparent that Cj in <33at is a more conve- 
nient variable than Cj in <34t because the former only depends only 
on the temperature whereas the latter depends on both the temper- 
ature (pressure) and number density. The results will only depend 
on the (physical) sound speed of one of the two species. 



4.3 Wave equations 



In the non-relativistic limit and using the definitions in the previous 
section, the evolution equations for the expansion and the auxiliary 
variables W and 'f are given by 

1 . 
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2Kp„,[A + m5] , 
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(35a) 
(35b) 

(35c) 



where in this limi t the expansion evolut ion <35at does reduce to 
equation (24b) in ISetschart et alJ J2004h for the gravitational at- 
traction of pressure-less dust. 

Since we want to study wave-like harmonic perturbations of 
the density and temperature we take the second time derivative of 
<29aH29d1 and use the evolution equations <35> to eliminate ^, Y 
and W from the second order differential equations for d. A, 7 and 

r 



A(T) = -[#(t) + ^(T)l, 
5(t) = -'f(T), 

r(T) = -|[#(T) + ar(T) + ^(t)l 



fi-c) 



[a#(T) + y(T) + aC(T)] 



(36a) 
(36b) 
(36c) 
(36d) 



Equivalently, we can simultaneously solve the combined set of 
first order ordinary differential equations <29> and <35L The solu- 
tions for initial conditions and limiting cases of interest are studied 
in the next section. 



5 SOLUTIONS 

In studying inhomogeneities as perturbations on a FRW Universe 
filled with an unmagnetized cold dusty plasma we expect the per- 
turbations to behave as eigenmodes of such a plasma altered by 
the background curvature. We will therefore study the longitudi- 
nal Langmuir mode, ion-acoustic oscillations and the gravitational 
Jeans instability. 
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5.1 General discussion 

By examining <35> we can make the following observations: 

• From <35at we find that an initial perturbation in either the 
total density or the charge density causes a (gravitational) pertur- 
bation of the expansion, unless the plasma consists of electrons and 
positrons with m = (in which case only a total density perturba- 
tion couples to the expansion); 

• A perturbation in the charge density corresponds to a longitu- 
dinal electric field perturbation (such that VgE" ^ 0) which enters 
the equations of motion as a Lorentz force term and excites per- 
turbations in the current density and (unless the plasma is e*) the 
total velocity variables. This is expressed in <35bH35ct in terms of 
the plasma frequency; 

• When the electric field perturbation is unimportant, only the 
wavelength dependent pressure gradient terms contribute in J3^}- 
I35ct . When the sound velocities of the two species are equal in the 
background, a perturbation in the total density only drives a per- 
turbation in the total velocity, while a perturbation in the charge 
density only excites a current density perturbation. If the plasma 
also were e^ => m = then <35bl435c1 completely de-couple and 
we have ^ = W{A,r) and "/ = 'f{S,Y)- However, in the general 
case that we want to investigate in this paper where the sound ve- 
locities in the background are allowed to differ, a perturbation in the 
charge density will excite a total velocity perturbation and more in- 
terestingly a total density perturbation will excite a current density 
from <35c> . 

The excited modes in C,, W and Y derived from <35> feed into <36> . 
or equivalently <29> and we can make some remarks on how they 
drive perturbations, particularly in the temperature variables. 

• The behaviour of the temperature variables can be found most 
readily by integrating the second expressions in <29cH29d1 and set- 
ting the integration constants to match the appropriate initial con- 
ditions. For a one-temperature background we find that the total 
temperature is simply proportional to the total density, whereas a 
temperature difference only couples to a perturbation in the charge 
density. 

• It again is clear from J36cll36d1 that a temperature difference 
in the background allows for additional coupling between the equa- 
tions. In particular, a perturbation in the total velocity or the ex- 
pansion excites perturbations in the temperature difference. For in- 
stance, a growing mode in the expansion such as gravitational col- 
lapse can lead to a growing temperature difference. 

This concludes our general discussion of the allowed excitations 
and couplings of the scalar perturbation variables in a non-trivial 
background. We now proceed to study some analytical solutions of 
these differential equations in more detail. 



5.2 Evolution of the background 

In solving the evolution equations for the perturbations we also 
have to include the evolution of the coefficients in the background. 
In the approximation that we introduced in S ection l4!2l where the 
pressure is included in the sound velocity but not in the gravita- 
tional mass we can just use the well known evolution of a dust 
dominated Universe. The scale factor then evolves as 5(t) oc r i 
and the dimensionless expansion and gravitational mass as = 2/t 
from (13) and kMN/4 = Kjl/l = 2/(3-r2) from (T2j. The num- 
ber density, which is hidden in the plasma frequency, evolves as 



Since the equation of state <8at is restricted to the non- 
relativistic regime we will not discuss the radiation dominated Uni- 
verse with p = p 13 here, although the extension is straightforward 
and could be the subject of a future paper. 



5.3 Long wavelength Langmuir modes 

The differential equations in the previous section are most readily 
solved in the long wavelength limit or alternatively the dust limit 
where f = everywhere. In both cases the terms proportional to 
T£ /S vanish. The e xact solutions derived i n this section for A and 
5 agree with those in lBetschart et alj (2004J) (where 5 is expressed 
in y = dn/N) but we can now extend these results to include the 
possible thermal effects. 

In this limit the wave equations for the charge density and total 
number density are 



8(t) 



A(t) 



4 S{t) 

"3~r~ 

4A(t) 
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2A(T) 



= 0, 
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(37a) 



(37b) 



The first mode we consider is excited by an initial perturbation in 
the charge density 5(1) = ^) when all the other perturbations are 
zero initially. From <29> we then find that 

A(l) = 5(l)=r(l) = y(I)=0. 

The solution of <37at subject to these initial conditions is 

5(t) = 5oT"5 [cos(a)lm:) + gLsin(coln-r)] , (38) 

where (O^ = O)^ — 1/36 ~ col. This oscillatory solution agrees with 



Bets chart et alJ <2004) . It is a non-propagating oscillation at the 
plasma frequency and is similar to the longitudinal Langmuir mode. 
The solutions are not simple plane waves but a more complicated 
harmonic oscillation with a decaying (t^''^) envelope due to the 
background expansion and corresponding evolution of the plasma 
frequency. 

On inserting the solution J37at in <37b> we can solve for A(t) 
to find 



with 



A(T)=m5o g(t) 



g(^) 



m 

* 



2,32 

-T ' + --r3 
5 5 



(39a) 



(39b) 



The solution <39> for the total density perturbation is a superposi- 
tion of the oscillatory Langmuir mode and the usual growing and 
decaying (power-law) modes of the standard gravitational insta- 
bility picture. The growing mode always dominates (since T > 1) 
leading to gravitational collapse on all length-scales. This is a well 
known result for dust since there is no pressure to act as a restoring 
force. 

The temperature perturbations can be found by solving the 
wave equations <36cH36dV but even simpler is to integrate the right- 
most expressions in J29cH29d1 taking care of the initial conditions. 



N/No = co'/co^ 



^lo 



r-2 



from n7bll23l . 



r(T) = 1^) 
r(T) = 1^) 



(a-"2)-T (a-mg(T)) 

(I-raa)^-— -(I-mag(T)) 



(40a) 
(40b) 
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From <40t we find that in an equal temperature plasma the rela- 
tive temperature perturbation only fluctuates with the oscillation 
in the charge density, but by allowing for a temperature differ- 
ence in the background (|o!| > 0) the temperature difference grows 
as a power-law during gravitational collapse. Note that from <33> 
(a — m) oc Sc, , so the total temperature perturbation only couples 
to the charge density oscillation when the two fluids have unequal 
sound velocities in the background (by comparison ( 1 — ma) oc Cj 
iniiObl). 



5.4 Gravitationally driven temperature difference 

We now study <37t in the same long wavelength or dust regime 
but start with the more common initial condition where the total 
density rather than the charge density is perturbed (so A(l) = Aq 
and all the other perturbations vanish at f = fo). With these initial 
conditions no perturbation of the charge density is excited and the 
total density only experiences the continuous gravitational collapse 
which couples to the temperature variables through the expansion 
<35> . The solutions are given in terms of the gravitational mode 
<39blbv 

5(t) = 0, (41a) 

A(t) = Aog(T), (41b) 

y(T) = §aAo(g(T)-l), (41c) 

r(T) = §Ao(g(T)-l), (41d) 

and we find from <41ct that even in this case gravitational collapse 
increases an initial temperature dijference when we start with a 
perturbation in the total density such as those that are observed in 
the cosmic microwave background (CMB). 

5.5 Wavelengtii dependent perturbations 

In treating perturbations of all wavelengths, the system of differen- 
tial equations becomes substantially more complicated. We there- 
fore specify to perturbations on time-scales that are small compared 
to the Hubble time such that the scale-factor is approximately con- 
stant ^(t) ~ 5o => = (or equivalently, the wavelengths are small 
with respect to the Hubble scale). Consequently, the Kp,„/2 terms 
vanish because of <12> and we can neglect the evolution of the co- 
efficients (A', cOp and f) in the background. In this limit, the pertur- 
bations locally behave as plane waves and we can solve the system 
of equations algebraically and derive a dispersion relation. 
For the wave equations we find in terms of ho = ^/5(1) 



(o}5-^ 



jKr\{^ f. 



Scl 
— =r< 



and from <29> we simply have 



'XtP-C^ 



y — raF 
1 —ma 
T — my 



1 —ma 



F=|[a5 + A] , 
y=|[aA+5] . 
A dispersion relation follows readily with solutions: 



co'pS, 


(42a) 
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(42b) 




(43a) 




(43b) 
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"p \ ^p 

If the ion-acoustic frequency is small compared to the total plasma 



frequency, we can expand the square-root term in k^ScJ cot, <C 1 
to find 



col ~ ^[c2 + |5c2|]=^max(c,,i,c,,2)^ 



(44a) 



CO 



2 ~ kj-Ac}-\5c^" ' -2 



'^0L>-.5 



co„ 



komm{cs^l,c,2) 



(44b) 



which are the ion-acoustic and Langmuir modes, respectively. 



5.6 Jeans criterion 

To obtain a Jeans length-scale for collapse, the gravitational inter- 
action is clearly essential. Therefore we cannot neglect icp,,, which 
means we cannot neglect and so we have to treat the full 
system of equations. The goal is to study the properties of the 
wave equation for the total density perturbations and find the criti- 
cal length-scale distinguishing collapse from oscillatory behaviour. 
From the non-relativistic limit of the expansion equation J35at and 
J36> one realises that rather than to study perturbations in the total 
number density A, a more physical choice is to study perturbations 
of the GI normalised density variable 



:A + m5 = 5- 



,V2p„ 



(45) 



where p,„ = ^ (MA' + 5m5n) to first order. Furthermore, we have 
to eliminate the temperature variables F and y from <35ll36t which 
is achieved by integrating <29t to find 

F — my = ^[{a — m)5 + {l—ma)A], 



y — mY = I [(a — m)A+(l — ma)5] 



(46a) 
(46b) 



Putting all of this together, we find 

p2 
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(l-m2) 



-m\Scj\)-^Kp„ 



Sc^lS. 



(47) 



We assume that the term on the right-hand-side of the equation is 
negligible for iri'^ I. Note also that the plasma frequency dependent 
terms have dropped out in the equation for S) which can be seen 
from OHEl or ^. 

Finding an exact solution of <47t is complicated by the fact 
that all the coefficients are themselves time dependent: 0(t), Cs{t), 
8cs{t), S(t) and Pm{'^)- However, a general result from differential 
calculus which in this particular case is called the Jeans criterion 
states that the solution for & changes from oscillatory to growing or 
decaying solutions when the term in square brackets changes sign 
(which corresponds to the perturbation frequency changing from 
real to imaginary). 

The limiting case is generally expressed as a Jeans wave- 
length Xj{t) — 2KS{t)/l which corresponds to the length-scale at 
which the plasma becomes unstable to gravitational collapse. For 
length-scales smaller than Xj the pressure prevents collapse. For an 
electron-ion plasma with m ~ 1 we find 

K 



^}- 



Gp, 



-{Q + \Scj\) 



(48) 



which corresponds to a well known result when Sc^ = (see for 
instance equation (220) in Ellis & van Fist 1 1999) ). In our more 
general treatment we find, given <33> . that in a two-temperature 
plasma, the Jeans wavelength is determine d by the species with the 
largest sound velocity. Note that, as in Ellis & van ElsJ il999h the 
Jeans length Xj is time-dependent. 
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The corresponding Jeans mass is given by 



M, 



AkX] 



;r] 2 4;rmax(cj.i,Cs2) 



m 



(49) 



The limiting cases discussed in Sections l5.3H5.4l are vahd for 
wavelengths much longer that the Jeans wavelength <48> . 



6 CONCLUSIONS 

In this paper we have extended the gauge invariant covariant pertur- 
bation theory on curved manifolds to include temperature effects in 
a multi-fluid that is not in complete thermal equilibrium. We have 
derived two closed sets of differential equations: one for gauge in- 
variant covariant vector quantities in Section|3]and one for corre- 
sponding scalar variables in Section |4| The system in Section |3| 
could be the subject of further investigations, but in this paper we 
have focused on studying the behaviour of the scalar perturbations. 

In the long wavelength limit for a pressure-less plasma we re- 
cover the results in Betschart et al. (2004) but we have generalized 
the treatment to a two-fluid where both species are perfect fluids 
with an ideal gas equation of state and find that the total temper- 
ature and relative temperature perturbations follow the total and 
charge densities, respectively. If the two fluids furthermore have 
slightly different temperatures in the background we find that they 
also couple to each other. Most interestingly, the temperature differ- 
ence grows during gravitational collapse. This result applies both to 
initial perturbations in the charge density or in the total density. We 
expect that the process of structure formation will be effected by 
this thermal evolution. 

To find propagating oscillatory modes we have included the 
wavelength dependent behaviour but neglected the evolution of the 
background which allows us to treat the perturbations as locally 
plane waves. In this approximation we have algebraically obtained 
a dispersion relation with eigen-frequency solutions corresponding 
to Langmuir and ion-acoustic modes. 

Finally, in Section l5!6l we have derived a modified Jeans wave- 
length that includes thermal effects and depends on the sound ve- 
locities of the two plasma species. In particular, the length-scale 
and mass for gravitational collapse are set by the species with the 
largest sound velocity. 

The theoretical framework presented in this paper could be 
used for a more detailed numerical analysis of the perturbations 
observed in the CMB, incorporating thermal effects in a gauge in- 
variant and covariant fashion. 



iBetschart et al.l y004J). The relations are linearized about a back- 
ground in which Wab = dah = "a = 0- X, Xi,, Xab Stand for any 
scalar function, vector or tensor and ^ can be any of those. 



{cvsVCgi,) = curlXoii, — ^0 curVC^b , 



VaV^X" = V^VaX" - I (0^ - 3^] WaX" , 

S cuxXiy'^X'') = V2(5 curUC") . 



(Ala) 
(Alb) 
(Ale) 
(Aid) 

(Ale) 
(Alf) 
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APPENDIX A: COMMUTATION RELATIONS 

In this Appendix we summarise the commut ation relations tha t are 
used throughout this paper and derived in iMaartensI il997» and 



